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Abstract 

The paper presents realization theory of discrete-time linear switched systems 
(abbreviated by DTLSSs). We present necessary and sufficient conditions for an 
input-output map to admit a discrete-time linear switched state-space realization. 
In addition, we present a characterization of minimality of discrete-time linear 
switched systems in terms of reachability and observability. Further, we prove that 
minimal realizations are unique up to isomorphism. We also discuss algorithms 
for converting a linear switched system to a minimal one and for constructing a 
state-space representation from input-output data. The paper uses the theory of 
rational formal power series in non-commutative variables. 

Keywords: hybrid systems, switched systems, realization theory, minimal 
realization. 

1 Introduction 

In this paper we develop realization theory of discrete-time linear switched systems 
(abbreviated by DTLSSs). DTLSSs are one of the simplest and best studied classes 
of hybrid systems, |f30| . A DTLSS is a discrete-time switched system, such that the 
continuous sub-system associated with each discrete state is linear The switching 
signal is viewed as an external input, and all linear systems live on the same input- 
output- and state-space. 

Realization theory. Realization theory is one of the central topics of system theory. 
For DTLSSs, the subject of realization theory is to answer the following questions. 

• When is it possible to construct a (preferably minimal) DTLSS state-space rep- 
resentation of the specified input/output behavior ? 

• How to characterize minimal DTLSSs which generate the specified input/output 
behavior ? 



*Univ Lille Nord de France, F-59000 Lille, France, and EMDouai, lA, F-59500 Douai, France, 
mihaly . petreczky @mines -douai . f r . 

tUniv Lille Nord de France, F-59000 Lille, France, and EMDouai, lA, F-59500 Douai, France, 
laurent -bakoQmines- douai . fr . 

♦Centrum Wiskunde en Informatica (CWI) P.O.Box 94079, 1090GB Amsterdam, The Netheriands 
J . H . van . Schuppen @cwi . nl 



1 



Motivation. While there is a substantial literature on linear switched systems, 
realization theory was addressed only for the continuous-time case 1201 [T9]| . The mo- 
tivation for devoting a separate paper to realization theory of discrete-time DTLSSs is 
the following. 

1 . Realization theory for DTLSSs is substantially different from realization theory 
for linear systems. 

2. Realization theory for DTLSSs is substantially different from the continuous- 
time case. More precisely, the realization problem both for continuous-time lin- 
ear switched systems and for DTLSSs can be transformed to the same realization 
problem for formal power series. The difference lies in the specific transforma- 
tion. 

3. Formulating reaUzation theory explicitly for discrete-time DTLSSs will be useful 
the identification of these systems. In fact, the results of this paper were akeady 
used in 11211 for analyzing identifiability of DTLSSs . 

Intuitively, the main difference between linear realization theory and that of linear 
switched systems is the following. For linear switched systems, the realization prob- 
lem is equivalent to the problem of representing a sequence of numbers (Markov- 
parameters) as products of several non-commuting matrices (pre- and post-multiplied 
by fixed matrices). For linear case, the corresponding problem involves not products 
of non-commuting matrices, but powers of one matrix. In addition, for linear switched 
systems we allow arbitrary non-zero initial state. The presence of a non-zero initial 
state means that the input response and initial-state response have to be decoupled. A 
similar approach was already described in [31] for linear systems. 

Contribution of tlie paper We prove that span-reachability and observability of 
DTLSSs is equivalent to minimality and that minimal realizations are isomorphic. We 
also show that any DTLSS can be transformed to a minimal one while preserving its 
input-output behavior, by presenting a minimization algorithm. In addition, we for- 
mulate the concept of Markov-parameters and Hankel-matrix for DTLSSs . We show 
that an input-output map can be realized by a DTLSS if and only if the Hankel-matrix 
is of finite rank. We also present a procedure for constructing a DTLSS state-space 
representation from the Hankel-matrix. Our main tool is the theory of rational formal 
power series |5 , 29|. 

Related work To the best of our knowledge, the results of this paper are new. The 
results on minimality of DTLSSs were akeady announced in llSTI . but no detailed proof 
was provided. The results on existence of a realization by a DTLSS were not previously 
published. 

The realization problem for hybrid systems was first formula ted in ifTTl . In ifTTl [35l 
the relationship between input-output equations and the state-space representations was 
studied. In It 1 8l |26l l22l realization theory for various classes of hybrid systems were 
developed. In particular, realization theory for continuous-time (bi)linear switched 
systems was developed in ll20l \T9 \. The approach of the present paper is similar to 
that of L201, however the details of the steps are different. There is a vast literature 
on topics related to realization theory, such as system identification, observability and 
reachability of hybrid systems, see lfT6ll6]|30]|2imi33]|Mll32l[T4ll27ll4ll8l[l31[Bin^ 
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Our main tool for developing realization theory of DTLSSs is the theory of rational 
formal power series. This theory was already used for realization theory of nonlinear 
and multi-dimensional systems, ||9] EH |29] El ■ State-affine systems from |29| include 
autonomous DTLSSs as a special case. Realization theory of state-affine systems is 
equivalent to that of rational formal power series. In this paper we reduce the realization 
problem for DTLSSs directly to that of rational formal power series. Hence, indirectly 
we also show that the realization problems for DTLSSs and state-affine systems are 
equivalent. One could probably reduce the realization problem for DTLSSs to that of 
state-affine systems directly, however it is unclear if such a reduction would be more 
advantageous. 

Outline f|2] presents a brief overview of realization theory of discrete-time linear 
systems. f|3] presents the formal definition of DTLSSs and it formulates the major 
system-theoretic concepts for this system class. ij4]- 52] states the main results of the 
paper fj6] contains the necessary background on the theory of rational formal power 
series. The proofs are presented in SjTjand AppendixlAl 

Notation Denote by N the set of natural numbers including 0. The notation de- 
scribed below is standard in automata theory, see IIToH TI. Consider a set X which will 
be called the alphabet. Denote by X* the set of finite sequences of elements of X. 
Finite sequences of elements of X are be referred to as strings or words over X. Each 
non-empty word w is of the form w — 0102 - ■■ ai^ for some 01,02, ■■■ ,0/^ e X. The el- 
ement Oi is called the ith letter ofw, for / = l,...,k and k is called the length w. We 
denote by e the empty sequence (word). The length of word w is denoted by |w|;note 
that |£| = 0. We denote by X+ the set of non-empty words, i.e. X+ =X*\ {e}. We 
denote by wv the concatenation of word w G X* with v e X*. We use the notation of 
ifTsl for matrices indexed by sets other than natural numbers . For each j = l,...,m, ej 
is the jth unit vector of M™, i.e. ej = {dij,... , 5„j), Sjj is the Kronecker symbol. 



2 Realization theory for linear systems 

In this section we present a brief review of realization theory of discrete-time linear 
systems, based on [311. Although the results of this section are not used in the paper, 
they help to get an intuition for the results on realization theory of DTLSSs . 

The input-output maps of interest are of the form y : (M'")+ RP. For each se- 
quence M = Mo ■ ■ ■ Mf , f > 0, y{u) is the output of the underlying system at time f , if inputs 
uq,... ,u, are fed. It is well-known that for y to be realizable by a linear system, it must 
be of the form 

r-l 

y{uo---Ut)^Kt + Y,Ht-j-iUj (1) 

j=o 

for some matrices K/i G K'', Hi^ G M.p^'", k = 0,l,2,..., and for any sequence of inputs 
uo,... ,u, G M'". Consider a discrete-time linear system 

^ I JCf+i — Axt + But where xq is fixed 
\yt= Cxt 
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where A, B and C are n x n, n x m and p x n real matrices and xq G M" is the initial 
state. Note that the initial state is xq, and xq need not be zero. The map y is said to be 
realized by E, if the output response of E to any input u equals y{u). This is the case if 
and only if y is of the form ([T]i, and Kt = CA'xq, Ht = CA'B, t >Q. We call E a minimal 
realization ofy, if it has the smallest state-space dimension among all the linear system 
realizations of y. 

Theorem 1 ( BID . Assume that E is a linear system realization ofy. Then E is a mini- 
mal realization ofy, if and only if it is weak- reachable and observable. Recall that E is 
weak-reachable if and only if {A, \xq B] ) is a reachable pair All minimal realizations 
ofy are isomorphic and any realization ofy can be transformed to a minimal one. 

The transformation to a minimal system can be carried out by first transforming the 
linear system to a weak-reachable one, and then to an observable one, |3r|. 

Next, we formulate conditions for existence of a linear system realization of y. 
To this end, we assume that y is of the form This assumption is necessary (but 
not sufficient) for existence of a realization. We call the matrices M, = \Kt fr], 
f > Markov parameters. This terminology is slightly different from the one used 
in II3TI . Note that y is completely determined by the Markov-parameters {M;}^q. In 
addition, note that we defined the Markov-parameters without assuming the existence 
of a linear system realization. In fact, we use the Markov -parameters for characterizing 
the existence of a linear system realization. More precisely, we define the infinite block 
Hankel-matrix Hy ofy as follows Hy = {Hi.j)'^j^i, Hi j = M,+ ,_2, i-e. the entries of Hy 
are formed by the entries of the Markov-parameters of y. 

Theorem 2 ( 1311 ). The map y can be realized by a linear system if and only if the rank 
ofHy is finite, //'rank Hy = n < +°°, then a minimal linear system realization E ofy can 
be constructed from the columns of Hy. In particular, this means that rank//,, equals 
the dimension of any minimal linear system which is a realization ofy. 

Procedure 1. The construction of Y, from the columns of Hy is as follows. Fix a finite 
basis in the column space of Hy. Then xo is formed by the coordinates of the first 
column ofHy in this basis, the rth column of the matrix B represents the coordinates of 
the r -\-\th column of Hy in this basis. The matrix C is the matrix (in the fixed basis) 
of the linear map which maps each column to the vector formed by its first p entries. 
Finally, A is the matrix (in the fixed basis) of the linear map which maps the jth column 
to the y + (m + \ )th column, i.e. it maps the block column (M,-|-y_2)^i to the block 
column (M,+y_i)~Li. 

3 Linear switched systems 

In this section we present the formal definition of DTLSSs along with a number of 
relevant system- theoretic concepts for DTLSSs . 

Definition 1. Recall from f21i that a discrete-time linear switched system ( abbreviated 
by DTLSS), is a discrete-time control system of the form 





(3) 
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Here Q = {I,. . . ,D} is the finite set of discrete modes, D is a positive integer. For each 
t ^N, qt G Q is the discrete mode, m, S M is the continuous input, yt G M.P is the output 
at time t. Moreover, G W^", € M"^'", Cq S W^" are the matrices of the linear 
system in mode q £ Q, andxQ is the initial continuous state. We will use 

{p,m,n,Q,{{A^,B^,C^) \ q E Q},Xq) 

as a short-hand notation for DTLSSs of the form (O. 

Throughout the section, E denotes a DTLSS of the form Q. The inputs ofT. are the 
continuous inputs {mjIJ^q and the switching signal {qt}'^^Q- The state of the system at 
time f is X,. Note that any switching signal is admissible. We use the following notation 
for the inputs of E. 

Notation 1 (Hybrid inputs). Denote = Q x W". 

We denote hy ^* (resp. '2'+) the set of all finite (resp. non-empty and finite) 
sequences of elements of . A sequence 

w^{qo,uo)---{q,,u,)e'i/^,t>0 (4) 

describes the scenario, when the discrete mode qi and the continuous input m,- are fed 
to E at time for i ~0,...,t. 

Definition 2 (State and output). Consider a state Xi„it £ R". For any w G of the 

form 0, denote by xz(xi„it,w) the state of E at time t +1, and denote by yz{xinitjw) the 
output o/E at time t, ifL is started from Xi„it and the inputs {m,};^q and the discrete 
modes {qiYi^Q are fed to the system. For notational purposes, we define xz{xi„it,£) = 

^init- 

That is, X£(x,„,,, w) is defined recursively as follows; x-£{xinit,£) = Xi„it, and if w = 
v{q,u) for some {q,u) g v G '^*, then 

XL{xinit,w) =AgXl,{xinit,v) +BqU. 

If w e and w = v{q,u), {q,u) G v e '^*, then 

yZ {Xinit , W) = CgXz {Xinit , v) . 

Definition 3 (Input-output map). The map yz'. ^ K'^ defined by "iw £ : 
3'l(w) — ^(xo, w), is called the input-output map ofL. 

That is, the input-output map of E maps each sequence w e ^2^+ to the output gen- 
erated by E under the hybrid input w, if started from the initial state xq. The definition 
above implies that the input-output behavior of a DTLSS can be formalized as a map 

f-.'W+^W. (5) 

The value /(w) for w of the form (|4]i represents the output of the underlying black -box 
system at time f, if the continuous inputs {m;}5=o switching sequence {qiYi=Q 

are fed to the system. This black-box system may or may not admit a description by a 
DTLSS. 

Next, we define when a general map / of the form (|5]l is adequately described by 
the DTLSS E, i.e. when E is a realization of /. 
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Definition 4 (Realization). The DTLSS T, is a realization of an input-output map f of 
the form (|5]l, if f equals the input-output map ofL, i.e. f =yz- 

The reachable set Reachijl) of E is the set of all states which can be reached from 
the initial state xq of E, i.e. 

Reach{I.) = {xz{xo,w) eR"\we^*} 

Definition 5 ((Span-)Reachability)). The DTLSS E is reachable, ifReach(Y,) — R", and 
E is span-reachable i/R" is the smallest vector space containing Reach{T.). 

Reachability implies span-reachability but in general they are not equivalent. 

Definition 6 (Observability). The DTLSS E is called observable if for any two states 

XuX2(zW of Z, 

(Vw e ^+ :y£(jci,w) =y£(jC2,M')) =^xi=X2 

That is, observability means that if we pick any two states of the system, then for 
some continuous input and switching signal, the resulting outputs will be different. 

Definition 7 (Dimension). The dimension ofY,, denoted by dimE, is the dimension n 
of its state-space. 

Note that the number of discrete states is fixed, and hence it is not included into 
the definition of dimension. The reason for this is the following. We are interested in 
realizations of input-output maps, which map continuous inputs and switching signals 
to continuous outputs. Hence, for all possible DTLSS realizations, the set of discrete 
modes is fixed. 

Definition 8 (Minimality). Let f be an input-output map. Then E is a minimal realiza- 
tion of/, ifL is a realization of f, and for any DTLSS E which is a realization of f, 
dimE < dimE. 

Definition 9 (DTLSS morphism). Consider a DTLSS Ei of the form ^ and a DTLSS 
E2 of the form 

E2 = (p,m,n",e,{(A«,B«,C«) | q G Q}A) 

Note that Ei and E2 have the same set of discrete modes. A matrix 5^ G R""^" is said 
to be a DTLSS morphism/rom Ei to E2, denoted by 5^ :Y^\ — ?> E2, if 

yxo = x^, and \fqeQ: .YAg, = YB,,, O^^Y = Q. 

The morphism Y is called surjective ( injective ) if ,Y is surjective ( injective ) as a lin- 
ear map. The morphism Y is said to be a DTLSS isomorphism, if it is an isomorphism 
as a linear map. 
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4 Main result on minimality 



Below we present the main results of the paper on minimality of DTLSSs. In addition, 
we present a minimization procedure and rank tests for checking minimality. In the 
sequel, E denotes a DTLSS of the form and f denotes an input-output map f : 

^+ -> RP. 

Theorem 3 (Minimality). 1. A DTLSS realization of f is minimal, if and only if it 
is span-reachable and observable. 

2. All minimal DTLSS realizations of f are isomorphic. 

3. Every DTLSS realization of f can be converted to a minimal DTLSS realization 
of f (see Procedure^below). 

The proof of Theorem|3]is presented in ij?] 

Remark 1. Note that can be minimal, while none of the linear subsystems is minimal, 
see Example\l\below. Since all minimal realizations are isomorphic, it then follows that 
such a DTLSS cannot be transformed to a one where at least one subsystem is minimal 
without loosing input-output behavior 

For analogous theorem for continuous-time linear switched systems see fSO", T^. 
Intuitively, the theorem says the following. First, a minimal DTLSS should not contain 
states which are not linear combination of the reachable ones (hence span-reachability). 
Second, a minimal DTLSS should not contain multiple states which exhibit the same 
input-output behavior (hence observability). Next, we present rank conditions for ob- 
servability and span-reachability. These conditions can be used to test minimality and 
to formulate Procedure |4] 

Notation 2. Let X be a finite set, ^ be a linear space, Aa : ^ ^ ,o £X be linear 
maps and let w € X*. The linear map A„ on ^ is defined as follows. Ifw = £, then Ag 
is the identity map, i.e A^x = x for all x G If w = <Ji<J2 - ■■ <Jk X*, (Ji, - ■ (Jk 'EX, 
k > 0, then 

A,;- ~ Afj^Aoj J • • • A(jj • (6) 

If ^ ~W for some n > 0, then A,v and each Aa, <J EX can be identified with annxn 
matrix. In this case A„ defines a product of matrices. 

We denote by Qf~" the set {w G g* | | w| < n} of all words w G g* of length at most 
« — 1. We denote by M„ the cardinality of Q^" and we fix an enumeration 

e<" = {yi,...,VM„}. 

We will use the notation defined above to define observability and reachability matrices 
for DTLSSs . 

Theorem 4. Span-Reachability. Define the span-reachability matrix R{T.) ofT. 

^(E)=[A,,B, A,,B, A,^„b] GK">^(lel»'+i)'^" w/iere 

B^[xo, Bu Bd] 
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Then E is span-reachable if and only i/rank ^(E) = n. 

Observability. Define the observability matrix (9(E) £ RPlGlMnX" o/E as follows. 



- CA,, - 








where C = 




-CAvM„- 




Cd_ 



Then E is observable if and only (/rank 0(E) = n. 

Informally, ,^(E) is formed by horizontal concatenation of blocks Ah-Z?,,, for all 
w G 2"^", q^Q, and (9(E) is the vertical concatenation of blocks QA,v, ^ G 2, w G 2^". 
Notice that if Q = {!}, then ^(E) is the controllability matrix of (Ai, [xq Bi]) and 
(9(E) is the observability matrix of (Ci,Ai). Hence, the linear system (Ai,Bi,Ci,xo) 
is weak-reachable (observable) if and only if it is span-reachable (observable), if inter- 
preted as a DTLSS. Hence, Theorem|3]implies Theorem[T] 

The result of Theorem |4] follow from |30|, the detailed proof can be found in Ap- 
pendix|A] Next, we formulate procedures for reachability, observability and minimality 
reduction of DTLSS s . 



Procedure 2 (Reachability reduction). Assume dim^(E) = n'' and choose a basis 
bi,. . . ,b„ ofW such that bi,. . . ,b„r span Im^(E). In the new basis, Ag,Bg,Cq, q G Q 
and XQ become as follows 



Aq = 



A'' A 

K 



0, 



: Cq \Cq , ^ ,Bq 







where A^ G 



G Mr Then E, = (p,m,«^e,{(AJ,B^,q) | q G 
Q},Xq) is span-reachable, and has the same input-output map as E. 



.B'qe\ 



Intuitively, E^ is obtained from E by restricting the dynamics and the output map of 
E to the space Im7?(E). 

Procedures (Observability reduction). Assume thatksxOil.) — n — n^' and let Z?} , . . . , bn 
be a basis in M" such that b„o+i, ■ ■ ■ span ker(9(E). In this new basis, Aq,Bq, Cq and 
xq can be rewritten as 



A" 







Aqi Aq 



c 

J ^q 



OJ ,Bq 



where A'q G 



r^" G 



C" G and x" G 



Then the DTLSS Eq = 



{p,m,n",Q,{{Aq,B'q,C'q) I q G 2}i-*^o) ''^ observable and its input-output map is the 
same as that ofL. IfL is span- reachable, then so is Eo. 



Intuitively, E^ is obtained from E by merging any two states xj, X2 of E, for which 
(9(E)xi — (9(E)x2. The latter is equivalent to yz{x\ , w) = yz{x2,w), Vw G 
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Procedure 4 (Minimization). First transform T, to a span-reachable DTLSS T,r and 
then transform to an observable DTLSS E,„ = {'^r)o- Then E,„ is a minimal realiza- 
tion of the input-output map ofL. 

The correctness of Procedures |2i3 | and|4]are proved in f|2l using the theory of formal 
power series. Note that the correctness of Procedure |3]and of Procedure |2] (in case of 
xq = 0) has already been shown by a direct proof in ll30l . 

Example 1. Let E = {p,m,n,Q,{{A^,B^,C^) \ q G Q],Xq) with Q = {1,2}, n 3, 
^o=[0 1 0]^ 



Ci = [1 0] 
C2 = [0 1] 



This system is observable, but it is not span- reachable. In order to see observability, 
notice that the sub-matrix \<3l {C\A\Y Cj] of 0{L) is of rank 2>. In order to see 
that E is not span-reachable, notice that if {x^y^zY ^ column of R{L), then z — Q- 
Hence dim7?(E) < 2. 

Using Procedure^ we can transform E to the minimal realization 

E„, = (p,m,n'",e,{(A- B'«,C-) | q £ Q],x'S) 
ofyz: Q = {1,2}, n™ = 2, x'^ = [l, O] ^ and 





"0 


1 


0" 




"0" 


Ai = 








1 


,Bi = 













1 









"0 


1 


0" 




"o" 







1 


1 


,B2 = 


1 










1 








A'l' = 





1 

1 0' 

1 



^Bi - 



,cr = [o, 1] 

,C?=[0, 0] 



Using 07V . it is easy to see that neither (A'",B'",C™,Xo ) nor (Aj j-B™ j^?'-^) '^'"^ 
minimal. 



5 Main results on existence of a realization 

We present the necessary and sufficient conditions for the existence of a DTLSS real- 
ization for an input-output map. In the sequel, / denotes a map of the form (|5]l. To 
this end, we need the notion of the Hankel-matrix and Markov-parameters of an input- 
output map. More precisely, we proceed as follows. First, we define the notion of 
Markov parameters of / and use them to define the Hankel-matrix of /. We then use 
the Hankel-matrix to formulate conditions for existence of a DTLSS realization of /. 
To this end, we need the following notation. 
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Notation 3. In the sequel, we identify any element w = (qo, uq) ■ ■ ■ {qt,ut) G with 
the pair of sequences (v, u), v G Q'^, u G (K™)+, v ~ qo - ■ -qt and u = uq - ■ ■ Uf. 

Notation 4. Consider the input-output map f. For each word v G of length |v| = 
t > define /v : {W")' W as 

Mu) = f{{v,u)). (7) 

Now we are ready to define the Markov-parameters of an input-output map. 

Definition 10 (Markov-parameters). Denote g*'* = {w G 2* | |w| > k}. Define the 
maps S{ : Q^ * W and sj : Q^ * -^W,j = l,...,mas follows; for any v G Q*, 

q,qo e Q, 

J- 

Soivq) =/v^(0,...,0) and 

f (8) 
Sj{qovq)=fqQvq{ej,0,...,0) -/^^^^(O, • • • ,0), 

with Cj G M™ is the vector with 1 as its jth entry and zero everywhere else. The collec- 
tion of maps {-S^IJ^o is called the Markov-parameters of f. 

f f 
The function can be viewed as the initial state -response and the functions , 

j = l,...,m can be viewed as input responses. The interpretation of 5q, sj will become 
more clear after we define the concept of a generalized convolution representation. 
Note that the values of the Markov-parameters can be obtained from the values of /, 
i.e. by means of input-output experiments. 

Notation 5 (Sub- word). Consider the sequence v = qo - • - qt d Q^, qo,---,qt ^Q,t>0. 
For each 7, A: G {0, . . . define the word Vj^^ € Q* as follows; if j > k, then Vj^^ = if 
j = k, then vjy = qj and if j < k, then vy|^ = qjqj+i ■ • • qh That is, Vj\i^ is the sub-word 
of V formed by the letters from the jth to the kth letter 

Definition 11 (Convolution representation). The input-output map f has a general- 
ized convolution representation (abbreviated as GCR), if for all w = (v,m) G 
V = qo---qt, u = uo---Ut, qo, - ■ ■ ,qt & Q, uo,...Ut G W", f{w) can be expressed via 
the Markov-parameters of f as follows. 

i-i 

f(w) = 4(vo|f-i?«) + L S^{qkVk+ii,-iqt)uk 

k=0 

where Sf{w)=[s{{w), s{{w), siiw)] eW'<'" for all w & Q*. 

Remarli 2. If f has a GCR, then the Markov-parameters off determine f uniquely. 

The motivation for introducing GCRs is that existence of a GCR is a necessary 
condition for realizabiUty by DTLSSs. More precisely, the following holds. 

Lemma 1. The map f is realized by the DTLSS Z if and only iff has a GCR and for 

all V G Q*, q,qo G Q, 

S'k (vq) = CaAyXQ and 

(9) 

(qovq) = CqAyBq^ej, j = l,...,m. 
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The proof of Lemma [T] can be found in Appendix |A] From Lemma [T] it follows 
that if / is realizable by a DTLSS, then the values of and sj, j = 1, . . . ,m can 
be expressed as products of matrices. Moreover, Sq corresponds to the part of the 
response which depends on the initial state, and {sj}^^^ encodes the response from 
the zero initial state. 

We can draw the following analogy with the linear case ^ Existence of a GCR is 
analogous to the requirement that the input-output map is of the form ([T]i. The Markov- 
parameter 5q(v^) corresponds to the vector K^y^, and the vector sj{qovq) corresponds 
to the jth column of the matrix //|,,|. Finally, if / can be realized by a DTLSS, then 
the Markov-parameters can be expressed as products of matrices (|9]l. This is analogous 
to the linear case, where K, = CA'xq and H, = CA'B holds for f > 0, if (A,fi,C,xo) 
is a realization of the input-output map. In fact,if Q = {I}, i.e. we are dealing with 



linear systems, then SQ{vq) — K^^^, S-^j{qQvq) is the jth column of and the GCR 

is the representation of the form dl}, and the right-hand sides of (|9]l becomes CA I' Ijcq, 
CAl'Ifie;, whereC = Ci,A =Ai,B = Bi. 

Next, we define the concept of a Hankel-matrix. Similarly to the linear case, the 
entries of the Hankel-matrix are formed by the Markov parameters. For the definition 
of the Hankel-matrix of /, we will use lexicographical ordering on the set of sequences 
Q*. 

Remark 3 (Lexicographic ordering). Recall that Q — {\,... We define a lexico- 
graphic ordering -< on Q* as follows. For any v,i £ Q*, v ^ s if either \v\ < \s\ or 
< |v| = \s\, V ^ s and for some I G {1,..., \s\}, v/ < i/ with the usual ordering of 
integers and v; ~ Sj for i — I, ... ,1 ~ I. Here v,- and Si denote the ith letter ofv and s re- 
spectively. Note that -< is a complete ordering and Q* = {v'l , V2, . . .} with v'l ^ V2 ^ 

Note that vi = e and for all i CzN, q Cz Q, v,- ~< Viq. 

In order to simplify the definition of a Hankel-matrix, we introduce the notion of a 
combined Markov-parameter 

Definition 12 (Combined Markov-parameters). A combined Markov-parameter M^^ (v) 
off indexed by the word v ^ Q* is the following pD x (Dm + 1) matrix 

'Slivl), Sf (Ivl), Sf{Dvl)' 



(v) = 



slivD), Sf{lvD), Sf{DvD) 



(10) 



Sm{w) 



where for any w £ g+, \w\ > 2, S^w) — s{{w), ^jCw), . 

Definition 13 (Hankel-matrix). Consider the lexicographic ordering -< of Q* from 
Remark^ Define the Hankel-matrix Hf of f as the following infinite matrix 



Hf = 



M'^(vivi), M'^(v2Vi), Mf{vkVi), 
Mf {viV2), M-f{v2V2), Mf{vkV2), 
M-f{viVi) Mf{v2V3), Mf{vkVj) 
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i.e. the pD x [niD +1) block of Hf in the block row i and block column j equals the 
combined Markov-parameter ivjVi) of f. The rank of Hf, denoted by rank Hf, is 
the dimension of the linear span of its columns. 

The Hankel-matrix of / can also be viewed as a matrix rows and columns of which 
are indexed by words from Q* . 

Remark 4 (Alternative definition of the Hankel-matrix). Notice that every row index 
<l ofHf can be identified with a tuple (v, /), / = 1 , . . . , pD and v G Q* as follows; 
V = Vr, i.e. V is the rth element of Q*, for some < r G N such that I = {r— I )Dp + /. 
In fact the identification above is a one-to-one mapping. 

Similarly, every column index < £ N can be identified with a pair [w,]) where 
w G Q*, j d Jf = {0} UQ X {I, . . . ,m}, where w — v,-, i.e. w is the rth element of Q* 
for some r G N such that k — (r — 1 ) {mD + 1 ) + i for some integer i — 1 , . . . , mD + 1, 
and if i = 1 then j — and if i = m{q —I) +z + I for some q £ Q and z ~ 1, . . . ,m, then 
j = {q,z). This identification is one-to-one. 

Using the identification of row and column indices outlined above, we can view Hf 
as a matrix, rows of which are indexed by (v, /), v G Q*, i — I,... ,pD, and columns of 
which are indexed by (w, J), w G Q*, j G Jf. The entry [///■] ^^^.^ of Hf indexed by 

row index (v, /) and column index [w^j) is the ith entry of the rth column of {wv), 
where r — I, if j — and r — m{q — I) + z + I if j = (q , z). In other words, 

[^/](,.,0.(w.(,,z))= [^^(^^^«')]/' 

where at = K + \ with K and I defined from i by the decomposition i — pK + I, K ~ 
0, 1 , . . . ,D — 1, / = 1 , . . . Here, [a]/ denotes the Ith entry of a vector a. 

It is not difficult to see that for Q — {\},Hfis the same as the Hankel-matrix defined 
in Sj2] The main result on realization theory of DTLSSs can be stated as follows. 

Theorem 5. The map f has a realization by a DTLSS if and only if f has a GCR and 
rank Hf < A minimal realization of f can be constructed from Hf (see Procedure 
15]) and any minimal DTLSS realization off has dimension rank Hf. 

Procedure 5. If rank Hf ~ n < +°o, then a DTLSS E/ of the form (O can be con- 
structed from Hf as follows. Choose a basis in the column space ofHf. 

In this basis, let xq be the coordinates of the first column of Hf. For each I = 
1, . . . ,m, the Ith column ofBq, q £ Q is formed by coordinates of the m{q — I) + 1 + Ith 
column ofHf. Let Cq, q € Q be the matrix of the linear map which maps every column 
to the vector formed by its rows indexed by p{q — \ ) + \ , p{q — I) + 2, . . . , pq. Define 
^q> q € Q as the matrix of the linear map which maps the rth column of the block 
column (M(vjV,))~Lj to the rth column of the block column {M{yjqvi))f^^, for each 
j = 1,2,..., and r — 1,2,..., (Dm + 1 ). 

Alternatively, using Remark^we can describe Hf as follows. The initial state xq 
is formed by the coordinates of the column of Hf indexed by (e,0). The Ith column 
of Bq, q (z Q is fondled by the coordinates of the column of Hf indexed by (e, {q,l)). 
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I = I, . . . ,m. The matrix Cq, q € Q is the matrix of the linear map which maps each 
column of Hf to the vector formed by its rows which are indexed by {£,p{q — 1) + 
1), . . . , {e,pq). Finally, Aq is the matrix of the map which maps each column indexed 
by (w,y) to the column indexed by {wq,j), w G Q*, j G Jf. 

Notice that for Q = {I}, Theorem |5] implies Theorem |2] and Procedure |5]reduces 
to Procedure [T] 

Example 2. Consider a SISO input-output map f such that for any v G Q^, \v\ = t, 

( 1 + L'fJi Uj ift > 1 and v = 2'-' 1 or 
/v(mi,...,m,) = < v = 2'-2ll, 

[ otherwise 

Hence, the Markov-parameters of f are as follows 

ift > 1 andv = 2'-h or v = 2'-^n 
otherwise 

ift >2 andv^ 2'-^l or v^2'-^ll 
otherwise 

It is easy to check that E from Example Q] satisfies (|9]l from Lemma |7] hence E is a 
realization of f. 

Consider the Hankel-matrix Hf of f. It is easy to see that the set of columns ofHf 
contains two elements: b\ and b2. The entries of bi equal I, if indexed by (v, 1) with 
\v\ > and v = 2l^'l or v = 21^''^' 1 and are zero otherwise. The only non-zero entry of 
b2 is 1 and it is indexed by (e, 1). Applying Procedure\5\to our example, and taking 
(^l,i>2) as a basis oflmHf, we obtain a DTLSS of the form (O which coincides with 
T.mfrom ExampleUl 

Indeed, since the column of Hf indexed by (e, (1, 1)) is zero, and the column in- 
dexed by (e,0) and (e, (2, 1)) is b\, we get Bi = Q, B2 — X(, = { \ Since the entries 
of any column indexed by (e,2) are zero, we get C2 = 0. Since the entries of bi and 
b2 indexed by (e, 1) are 1, we get Ci — (1, 1)^. Note that if the column of Hf indexed 
by {w,j) equals bi, then the column indexed by (wl,y') equals b2, the column indexed 
by (w2,y') equals b\ +b2. If the column indexed by {w,j) equals b2, then the column 
indexed by {wl,j) and (w2,j) are both zero. Hence, if Ai and A2 are viewed as linear 
maps on ImHf, then Aibi ~ b2, Aib2 — 0, = 0, A2b\ = b\-\-b2. In other words, 
the matrices Ai andA2 are precisely the same as the matrices A"^ and A2 from Example 

m 

Note that once the Markov-parameters are defined, the definition of Hankel-matrix 
presented above coincides with that of the continuous-time case. As a consequence, 
we can repeat the realization algorithm described in (25. Algorithm 1] for DTLSSs . 
Moreover, [,25 ^ Theorem 4] holds for DTLSSs . For the sake of completeness, below 
we state the reaUzation algorithm and its correctness explicitly for DTLSSs . 

Definition 14 {Hf im sub-matrices of Hf). For L,M G N define the integers II = 
N{L)pD and Jm = N(M)(;7iD+ 1). Denote by Hfi,M the following upper-left II x Jm 
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sub-matrix ofHf, 



Mf (vivi), 
Mf{viV2), 



M-f (V2V2), 



Algorithm 1 

Inputs: Hankel-matrix//^jvjv+i- 
Output: DTLSS Ea, 



1: Compute the decomposition Hf,N,N+i = OR such that O G R^'^xm r g i^«x/a,+i 

and rank R = rank O = n. 
2: Consider the decomposition 



R 



■'''N(JV+1) 



such that Cv, G M"x(0'"+i). Define R,R^ e W-''^, q e Q as follows 



R 



R„ 



c c 



Cv 



3: Construct Eat of the form ^ such that 

[xo,Bi,...,Bo] = 

the first mZ) + 1 columns of R 

[Cf , Cj , . . . , C|;] ^ = the first pD rows of O 

V?ee:A^ = R^R+, 

where R^ is the Moore-Penrose pseudoinverse of R. 
4: Return Ln 



(11) 

(12) 

(13) 



Remark 5 (Implementation). One way to compute the factorization Hf,N.N+i = OR 
is as follows. If Hf N N^i = t/EV is the SVD decomposition of Hf ^j^^i, then define 
= f/Ei/2 fl«t/R = Ei/2y. 

Theorem 6. //'rank Hf j^j^ = rank i/^, f/zen f/ze algorithm returns a minimal realization 
of f. The condition rank Hff^^f^ = rank Hj holds for a given N, if there exists an DTLSS 
realization E off such that dimE < 1. 

The proof of Theorem|6]can be found in ijTl 

Remark 6 (Computation of Hf^N^N)- Note that Hf^N^N can be computed from the re- 
sponses of f. However, in principle, the computation ofHj fj fj requires an exponential 
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number of input/output experiments involving different switching sequences. This is 
clearly not very practical. It would be more practical to build Hf^isf^ff based on the re- 
sponse of f to a single switching sequence. Preliminary results on the latter approach 
can be found in l[23]l . A detailed discussion of this approach goes beyond the scope of 
this paper 

6 Formal Power Series 

In this section we present an overview of the necessary results on formal power series. 
The material of the section is an extension of the classical theory of |5 29J, for the 
proofs of the results of this section see ifTS] l20l . 

Let X be a finite set, which we refer to as the alphabet. A formal power series S 
with coefficients in M'' is a map 

S:X* ^ W' 

We denote by M.'^ ^ X* ^ the set of all such maps. Let J be an arbitrary (possibly 
infinite) set. A family of formal power series in ^ X* indexed by J, abbreviated 
as FFS is a collection 

^ = {Sj e r' < X* >| ./■ e /}. (14) 

In the sequel ^ denotes a FFS of the form d7?l ). Notice that we do not require Sj, j €J 
to be all distinct , i.e. S[ — Sj for some indices J, / GJ,j^lis allowed. 

Let J be an arbitrary set and let li > 0. A d-J rational representation over the 
alphabet X is a tuple 

/? = (^,{A^}„ex,B,C) (15) 

where ^ is a finite-dimensional vector space over M, for each a G X, Aa : 5C ^ 
is a linear map, C : ^ — > M'^' is a linear map, and B = {By G J?^ | y G 7} is a family 
of elements of X indexed by J . If d and J are clear from the context we will refer 
to R simply as a rational representation. We call the state-space , A^, (7 G X the 
state-transition maps, and C the readout map of R. The family B is called the family of 
initial states ofR. The dimension dim ^ of the state-space is called the dimension of 
R and it is denoted by dim/?. If ^ = R", then we identify the linear maps A^, (J G X 
and C with their matrix representations in the standard Euclidean bases, and we call 
them the state-transition matrices and the readout matrix respectively. 

The d — J representation R from ( fTSl l is said to be a representation of^, if 

Vj G /, Vw G X* : Sj (w) = CA,,Bj, (16) 

where Notation |2] has been used. We say that the family *P is rational, if there exists 
a d-J representation R such that R is a representation of W. A representation R,„i„ of 
is called minimal if for each representation R of ^, dim7?„„„ < dimR. Define the 
subspaces 

Wr = Spm{A,yBj e ,Sr \weX*,\w\<nJ eJ} (17) 
Or = n kerCAv.. (18) 

wSX* ,|w| <n 
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We will say that the representation R is reachable if dimWR = AxmR, and we will say 
i\\?iiRKobservableifOR^{Q).hQiR^{^,{Aa}aex,B,C),R^{i^,{Aa}aex,B,C) 
be two d — J rational representations. A linear map .5^ : ^ — >^ ^ is called a represen- 
tation morphism, and is denoted by ^ : 7? — if 

yAa^Aa^ya ex, s^Bj = Bjyj ej, c = c.y (i9) 

If is bijective, then it is called a representation isomorphism. If is an isomor- 
phism, then R and R are representations of the same FFS , and R is observable (reach- 
able) if and only if R is observable (reachable). 

Remark 7. Let R be a representation of^' of the form ([15i , and consider a linear iso- 
morphism S^-.SC^ R", n = dim/?. Then YR = (R",{^Aa^"'}aGX,^fi,CJ5^"'), 
where S^B = {yBj G R" | y £ 7} is a representation of^' and it is isomorphic to R. 
The representation S^R is defined on an Euclidean space and its state-transition and 
readout maps can be viewed as matrices. 

Definition 15 (Hankel-matrix). Define the Hankel-matrix Hm of^ as the infinite ma- 
trix, the rows of which are indexed by pairs {v,i) where v G X*, i — \ ,... ,d, and the 
columns of which are indexed by [w,]) where w £ X*, j G J. The entry [^'j'](iy)^(H,j) of 
Hvfi indexed with the row index (v, /) and the column index (w, j) is defined as 

[^>i'](v,-)(v,./) = [^./(v^^)], (20) 

where [Sj{wv)\. denotes the ith entry of the vector Sj{wv) G R^'. The rank of Hm is 
the dimension of the linear space spanned by the columns of H\^, and it is denoted by 
rank //>j/. 

Tlieorem 7 (Existence and minimality, lfT8ll20l ). 1. The family I' is rational, if and 
only i/rank H^i < +°o. 

2. Tjfrank Hvy < +°°, then a minimal representation R of^ can be constructed from 
H\^, see Procedure^ 

3. Assume thatRmtn is a representation of^. Then Rmin is a minimal representation 
of^', if and only if Rmin is reachable and observable. If Rmin is minimal, then 
rank //ip = dim/?m,„. 

4. All minimal representations of^* are isomorphic. 

5. Any representation Rof^* can be transformed to a minimal representation Rmin 
of^, see Procedure^ 

We conclude by presenting procedures for reachability and observability reduction, 
minimization of representations and construction of a representation from the Hankel- 
matrix. In the sequel, /? is a representation of ^ and R is of the form ( fTsl l. 

Procedure 6 (Repr from Hankel-matiix, ([T8]|20|). //rankZ/ip < +°°, then 

R^ = {lTaH^,{Aa}cex,B,C) 
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is a representation ofW. Here, for each G & X, is the linear map which maps every 
column ofH^i indexed by (w, j) to the column indexed by {wO ,j). The initial states are 
B = {Bj I j G J}, where Bj is the column ofH^i indexed by {£,j), j G J- Finally, C is a 
linear map which maps every column ofH^i to the vector fanned by those rows of this 
columns which are indexed by (£, 1), . . . , (e,c/). Recall that R'' is set of coefficients of 
the formal power series Sj of^, j G J, i.e. Sj : X* M''. 

Procedure 7 (Reachability Reduction). Assume R is a representation of and it is 
of the form ( I15I ). Recall the definition of the reachable subspace Wr of R from f l77l ). 
Define the representation Rr — {Wr, {A'fj}aex,B'' jC^), where for each (J & X, is the 
restriction ofA(, to Wr, B' — {B j ^ \ j G J} = B, and C is the restriction ofC to 
Wr. Then R, is a reachable representation of^'. 

Procedure 8 (Observability Reduction). Assume R is a representation of^ and it is of 
the form (115b . Recall from (\18]l the definition of the observability subspace Or. Define 
the representation Rq — /OR^,{Aa}aex,B,C). Here ^ /Or is the quotient space 
of 3(y with respect to Or. Denote by [x], x <G 56' the equivalence class of all those 
y G !X such that x — yG Or. Then Ao-H = \Acx\, a G X, C\x\ = Cx for all x G 
and B = {Bj G St^ /Or | 7 G 7} is such that Bj = [Bj], j G J. Then Rq is an observable 
representation ofai' and ifR is reachable, then so is Ro. 

Procedure 9 (Minimization). A representation Rof^* can be converted to a minimal 
representation as follows. Use Procedure [7| to obtain a reachable representation R^. 
Apply Procedure\8\to Rr and obtain the observable representation Rmin — {Rr)o. Then 
Rmin is a minimal representation of^*. 

If J is finite, then Procedures|6]|2]|8] and|9]can be implemented, see ifTsll . 

More precisely, we can formulate a realization algorithm for rational representa- 
tions, f24l. Below we present slight extension of the results of fSSl |29l IT2I on real- 
ization algorithms for formal power series. The proofs of the results can be found in 
lUSllSll. We introduce the following notation. Let K,M efi. 

Im = { (v, /) I V G v| < M, / = 1 , . . . , p] 
iK = {{^vJ)\i^J,w^X\\w\<K} 

Intuitively, the elements of 1m (resp. Ja:) are those row (column) indices of i/ip, the 
X*-valued component of which is of length at most M (resp. K). 

Definition 16. Define the matrix H^ m.k os the matrix, rows of which are indexed by 
the elements of 1m, columns of which are indexed by the elements of Ik, and its entry 
(^>P,m,a:)(v,!),(wj) indexed by the row index {v,i) G 1m ond the column index {w,j) G 
Ik is defined as (77>p.m,a:)(v,(), (»,,/) = (^1')(v,i),(h'j) = {Sj{wv))i. The rank ofH^^M,K, 
denoted by rank Hx^m.k, is the dimension of the linear space spanned by its columns. 

That is, Hm M.K is the sub-matrix of //>j/ formed by the intersection of the columns 
indexed by the elements of Ik and of the rows indexed by the elements of 1m ■ If 7 is 
finite, then Hyfi M,K is ci finite matrix. 
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Theorem 8 (Realization algorithm, f24]). /frank H^in.n = rank H<ii, then the represen- 
tation Rn, to be defined below, is a minimal representation of^. Tfranki/vp < 1, 
then rank — rank holds. The representation is of the form ( I-/5D , with 

the state-space 2j = Imf/\j<jv,A'+l. c>-nd such that if we denote by C^j the column of 
ImHxfi indexed by {w,j) G In, then 

y{wj) e In ■ Aa{Cn.j) = Cn-aJ 

V(w,;)eJ;v:C(C.,j)= [C„,y((e,l)), C^,j{{E,p))f 
V; G y : Bj = C^,,- 

Here Cwj{{£,i)) is the entry of the column C^.j indexed by (e,/). i.e. it equals 
iH^,N,N){e,i),{w,jy i = 1, ■ ■ ■ 



7 Proof of the main results 

The proof of the results on realization theory relies on the relationship between formal 
power series representations and DTLSSs state-space representations. This relation- 
ship is completely analogous to the one for linear switched systems in continuous time, 

121 mi. 

Consider an input-output map / and assume that / has a GCR. Below we define 
the FFS ^' f associated with f. We also define the representation R^ associated with 
a DTLSS E and a DTLSS associated with a rational representation R. These no- 
tions allow us to relate FFS and input-output maps and to relate DTLSS with rational 
representations. In turn, these correspondences enable us to translate the realization 
problem for DTLSS to the problem of rationality of FFS. 

We first define the FFS associated with /. To this end, recall the definition (HJ of 
the Markov-parameters of /. 

Definition 17 (FFS associated with /). For each q eQ, each index j = l,...,m, define 
the formal power series S^,y,So G 2* ^ as follows; for each word w G Q*, 

discrete mode q (z Q and index j = 1 , . . . , m. 



h,.j){^)^[{S■fJ{q^^l))\ (5f(^w2))^ {sf^{qwD)f^ 

{si{w\)Y, {si{w2)Y, {si{wD)Y^ 



(22) 



So(h') 

Let Jf = {0} U{{q,l) I q G Q,l = 1, . . . ,m} and define the FFS associated with / by 

^/ = {§; G < e* > I G 7/} . (23) 

Notice that the values of S(^j)(vv) and Sq{w) are obtained by stacking up the 

Markov-parameters of sj{qwi) and 5'q(w/) respectively, for i — 1,...,D. Next, we 
define the representation associated with E. 

Definition 18. Assume that E is of the form Define the representation 7?£ associated 
with 'Las a p\Q\~Jf representation of the form ( fTSl l, where Jf = {0} U 2 x {1, . . . ,m} 
and the following holds. 
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• The alphabet X ofRz is the set of discrete modes Q, and d ^ p\Q\. 

• The state-space 3^ ofRz is the same as that ofL, i.e. ~ M". For each q G Q, 
the state-transition matrix A ofRz is identical to the matrix A^j ofL. 

• The p\Q\ X n readout matrix C is obtained by vertically "stacking up " the matri- 
ces Ci, . . . ,0), i.e. 

c=[cf, cj, clY 

• B = {Bj G ^ \ j G Jf}, where Bq = xq andBi^^ ;) is the Ith column of the matrix 
B, ofZ. 

The intuition behind the definition of Ry_ is that we would like R^ to be a represen- 
tation of if and only if (|24] | holds. Then the Aq matrices of the representation Rz 
should coincide with the Aq matrices of E. The initial states of Rz should be formed by 
the vector Bq (in order to generate So), and B^ej (in order to generate ^{q.j))- Finally, 
the readout map C should be formed by "stacking up" the matrices Q. Next, we define 
a DTLSS 'Lr based on a representation R. 

Definition 19. Consider a p\Q\~ Jf representation R of the form ( I15l l, over the alpha- 
bet X = Q with d — p\Q\. If ^ = W does not hold, then replace R with the isomorphic 
copy .yR defined in Remark\7\whose state-space is R". In the rest of the construction, 
we assume that ^ =M" for n — dim ^ holds and that Aq, q G Q are n x n matrices, 
andC is a p\Q\ x n matrix. Define the DTLSS associated with R as follows. Let Y-r 
be of the form di]) such that 

• for q €z Q, the matrix Aq ofY-R is identical to the state-transition matrix Aq ofR. 

• For each q (z Q, the matrix Cq is formed by the rows {q ^ l)p + l,{q — l)p + 
2,...,qpofC,i.e. ^ 

C = [C[ , , • • • , Cq] . 

• For each q E Q, Bq — [^(,^,1) , • • • 5(^,m)] • The initial state xq ofLR is defined 
as XQ ~ Bq. 

The intuition behind the definition of T.r is the following. We would like T.r to be 
such that if we apply Definition [18] to it, then the resulting representation Rz^ should 
be close to R. 

The relationship between the various concepts introduced above is as follows. 
Tlieorem 9. 1. The Hankel-matrix H^i^ equals the Hankel-matrix Hj of f. 

2. The representations R and Rz^ are isomorphic, and T-r^ = E. 

3. The DTLSS E is a realization of the input-output map f if and only if the associ- 
ated representation Rz is a representation of^'f. 

4. The representation R is a representation of ^' j if and only if the associated 
DTLSS T-R is a realization off. 
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5. The DTLSS E is a minimal realization of the input-output map f if and only if 
the associated representation Rz is a minimal representation of^* f. 

6. The representation R is a minimal representation of^'f if and only if the associ- 
ated DTLSS Y-R is a minimal realization of f. 

7. The DTLSS E is span-reachable (observable) if and only if the associated repre- 
sentation R-£ is reachable (resp. observable). 

8. The representation R is reachable (observable) if and only if the associated 
DTLSS "Lr is span-reachable (resp. observable). 

9. Assume that Y.\ and E2 two DTLSSs with the state-spaces R" and M"" re- 
spectively. A matrix ,5^ G R""^" is a DTLSS morphism J?^ : Ei — > E2 if and only 
if ^ : — > Rz2 is a representation morphism, if is interpreted as a linear 
map. 

The statements of Theorem|9] above are summarized in Table [T] 

Proof of^ Proof of Part [H Straightforward. 
Proof of Part|2j Straightfoi-ward. 

Proof of Part |3] and Part|4j The proof is analogous to the proof of Theorem 10 
from ||20| . First, note that if 7? is a representation of then R satisfies the assumptions 
of Definition [191 Since R is isomorphic to Rz^, Part |4] follows from Part|3] Part|3] 
follows by noticing that E is a realization of /, if and only if for all <7o G 6' ./ = 1 1 • • • : 

S(./oJ)W = [cr, Cj, Ca^A,,B,„e,and 
§o(w)=[Cf, Cj, C^] A„XQ. 

The above statement follows from Lemma [T] by taking into account the definition of 
§0 and But ( l24b is equivalent to /?£ being a representation of ^ f. Indeed, the 

matrix \C\ Cj , ■ ■ • C^] ^ in the right-hand side of (l24l i equals the readout matrix 
C of and the vectors B^^^Cj and jcq coincide with the initial states and Bq of 

Rz- Hence, (|24li in fact says that ^j(w) = CA„Bj for all w e Q*, j G 7/, i.e. that Rz is 
a representation of 

Proof of Part |5] and Part |6l Follows from Part |3] and Part 5] by noticing that 
dimE — dim7?£ and dim/? = AivciLr. 

Proof of Part |7] and |8l Since R^^ is isomorphic to R, it is enough to prove Part 
I2I To that end it is enough to show that Wr^ = Im/?(E) and Or^ = ker(9(E), i.e. the 
image of the reachability matrix of E equals the space Wr^ of Rz, and the kernel of the 
observability matrix of E equals Or^. 

Assume that Rz is of the form ([BJ, with ^ = R",d = p\Q\ andX = Q. To see that 
Im/?(E) = Wflj., notice that Im^(E) is the linear span of the columns of matrices An,Bq 
and vectors A^xq, q G Q,w € Q*, \w\ < n. But the initial states B of Rz consists of the 
columns of the matrices B^, q G Q, and of the vector jcq- Hence, Im^(E) is spanned by 
vectors A^Bj, j e Jf and hence it equals Wr^. 
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Similarly, the kernel of 0(E) equals the intersection of kerC^A,^, q £ Q,w £ Q*, 
\w\ < n. It is easy to see that H^ggkerCgAn. = kerCA,y, hence, ker(9(E) is the intersec- 
tion of all spaces kerCAn-, w G Q*, \w\ < n. But the latter intersection equals Or^. 

Proof of Part 121 The proof is analogous to the proof of Lemma 10 of ll20l . Since 
the state-spaces of R^^ and Ei are the same, and the state-spaces of R^^ and E2 are the 
same, ,^ can indeed be viewed both as a potential representation morphism from 
to R-£^ and as a potential DTLSS morphism from Ei to E2. Then it is enough to prove 
that satisfies (O with R = R^i and R = Ri,^ if and only if 5^ satisfies Definition^ 
The latter proof is routine. Indeed, assume that Ei is of the form ^ and that E2 is of 
the form 

E2 = («',e,{(A;,,B;,,c;,) i<?ee},x;). 

Assume that/?!, is of the form (fTSll and Ry-. = (M"' , {A^}<,eg,B',c') where b' = {b'j 

j G Jf}. Note that the matrices Ag and Ag of R^^, respectively coincide with the 
corresponding matrices of Ej and E2. Then isa DTLSS morphism if and only if 

(V^ e e : S^Ag - a[^S^, Cg = C'gy, yBg = Bg) 

and o5^xo = Xq . 

But yq ^ Q : Cq — Cg5^ \& equivalent to C = C ,y, since 

Similarly, yBg = b'^ is equivalent to: V/ = 1, . . .,m,yBqei — yBf^gj-^ =B'^ei =B'^^iy 

This, together with ^xq = x'q, implies that ^Bj — Bj for all j e Jf. 

Hence, we have established that y is a DTLSS morphism if and only if\/q^Q: 
yAg = A'^y, C = C'y, and yjeJf. yBj = B'y But this means that y-.R^^^ Rzi 
is a representation morphism. □ 

Proof Theorem\3l By Theorem|9] Part|5] E is a minimal DTLSS realization of / if and 
only if R^Rz is minimal. By Theorem]?] R is minimal if and only if R is reachable and 
observable. By Theorem |9] Part|7] the latter is equivalent to E being span-reachable 
and observable. Next, we show that minimal DTLSS realizations of / are isomorphic. 
Let E and E be two minimal DTLSS realizations of /. By Theorem|9] Part|5] Rz and R^ 
are minimal representations of . Then from Theorem|7]it follows that there exists a 
isomorphism y :R^-^ R^. From Part|9]of Theorem|9]is then follows that ^ : E — > E 
is an isomorphism. Finally, the correctness of Procedure 2]is shown in Remark[TO] □ 

Proof of Theorem\5\ Necessity 

Assume that E is a DTLSS which is a realization of /. Then by Lemma [l] / has a 
GCR. Moreover, by Theorem |9] Rz is a representation of 'P/, i.e. ^/ is rational. By 
Theorem|9] Part[T] and Theorem |7] the latter implies that rank Hf < +°o. 
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Sufficiency 

Assume that / has a GCR and rank Hf < +00, Then by Theorem|9] Part[T] and Theo- 
rem|2] is rational, i.e. it has a representation R. Then by Theorem|9]the DTLSS Y-r 
is a realization of /, i.e. / has a realization. 

Finally, the correctness of Procedure |5]follows from Remark [TT]below. □ 

Now we are ready to analyze Procedure 12 J3J l4]and 151 

Remark 8 (Correctness of Procedure |2). Procedure^is equivalent to the following 
procedure. Apply Procedure^to to obtain R^. Then T-y from Procedure\2\and 
are isomorphic. It then follows that is span-reachable, since Ry is reachable, and 
E,- and E have the same input-output map, since both R^ and Ry are representations of 

Remark 9 (Correctness of Procedure |3). Procedure\3\is equivalent to the following 
procedure. Apply Procedure |S] to Rz to obtain an observable representation Rg. It 
follows that Y-o from Procedure\3\and Er^ are isomorphic. Since R,, is observable, Eo 
is observable as well. IfL is span-reachable, then Rz is reachable. Hence, then Rg is 
reachable and thus Eo is span-reachable. Finally, both R^ and Rg are representations 
of ^'y^, from which it follows that the input-output maps ofY, andY,,, coincide. 

Remark 10 (Correctness of Procedure H)- Procedure^can be restated as follows. 
Apply Procedure |9] to Rz and denote the resulting minimal representation by R„,- It 
then follows that Ymfrom Procedure^is isomorphic to Er„,. Since by Theorem\9\T,R^^ 
is a minimal realization ofyz, then so is E^. 

Remark 11 (Correctness of Procedure |5). Procedure\5\can be reformulated as follows. 
Use Procedure^ to construct a minimal representation R of ^' f from Hf = Ht^^ . Then 
by Theorem^ Hr will be a minimal realization of f. It is easy to see that the DTLSS 
Hf from Procedure\5\is isomorphic to E^. 

We will continue with the proof of Theorem |6] 

Proof of Theorem^ The proof is almost the same as that of the continuous-time case, 
described in |25|. From Theorem|9]it follows that Hf K.L coincides with H-^^x.l, and 
hence, rank Hf ^ ^ — rank Hf is equivalent to rank H\fi^ ^ ^ = rank i/>p^. 

Assume now that rank Hfj^j^ — rank Hf. Then the representation R^ from Theo- 
rem[8]is well-defined and it is a minimal representation of Consider Algorithm[T] 
and the decomposition defined there. Then IrnHf ^N+i — ImO and there exists a left 
inverse 0+ G M"^''~ of O such that O+O = /„. 

Consider the linear map : lmHf^N,N+i where S^{x) = O^x for all x G 

ImHf N.N+i and recall that ///;v,w+i = HWf.N.N+i- It then follows that is a linear 
isomorphism, and its inverse is O. Moreover, the isomorphic copy 

S'Rn = (M",{J?^A,^-i},ee,{^(By) I e Jf},CS^-') 

of Rn is also a minimal representation of^'f. 

Consider now the DTLSS E.j^a^ associated with S^Ri\i. It is easy to see that the 
DTLSS E satisfies (II lH13l l and hence it coincides with the DTLSS Eat returned by 
Algorithm[T] Theorem|9]it follows then that Eat is a minimal realization of /. 
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Assume that there exists an DTLSS realization L of /, such that dimE < 1. 
Then by Theorem |5] rank Hf = rank H^ij < dimE < A' + 1. Hence, by Theorem [8] 
rank i/ip^ = rank //ij/^.jv.A'- D 

We conclude this section with the following remark. 

Remark 12 (Continuous-time case). If instead of a discrete-time system we consider 
a continuous-time system E, then the constructions ofR^ and T-r are exactly the same. 
The construction of^f differs only in the way the Markov-parameters sj{qQvq) and 

^0 {^4)' ^ G Q*' ItIo & Q, j — ^, ■ ■ ■ jin, are derived from the input-output map f. How- 
ever, S^iyq) = CqAyXQ and S'^-{qovq) = CqA^Bq^ej also holds for the continuous -time 
case, ifL is a realization of f. A detailed description of the continuous-time case can 
be found in l[2U\\T9]l. 



Realization of / 


Representation of ^/ 


j: = e„, <= 


Rz 


Er <= 


=> R = Rl, 


observable, span-reachable <S= 


^ observable, reachable 


minimal <^ 


^ minimal 


y, DTLSSmorphism <= 


y, representation morphism 



Table 1: Correspondence between DTLSSs and representations 



8 Conclusions 

We presented realization theory for discrete-time linear switched systems. The results 
and the proof techniques resemble the ones for continuous-time linear switched sys- 
tems presented in our previous work. 
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A Technical proofs 



Proof of Lemma\J] Consider the input-output map yz of E. By induction on t, it fol- 
lows that if w = (v, m) (1'^'^ ,v ~ qo - ■ ■ qt, u ^ uq - ■ ■ Ut,t >0,qo,...,qt EQ,uo,...,Ut € 
W", then 

t-i 

yziw) = Q,A,,o|,_,xo + CciAvj+y^t-^B^^Uj. (25) 

Consider the Markov-parameters S^{sq), S^iqasq), q,qo E Q, s E Q*, j = l,...,m, 
of yz- It then follows from (|25j and the definition of Markov-parameters that for all 
sEQ*, 

Si"" (sq) = QA,xo and Sf {qosq) = QA,B,,oey. (26) 

Notice that (l25T l - ( |26] | implies that yz has a generalized convolution representation. 

Assume that Lisa realization of /. Then yz = /. Then from (l25Tl-(l26]l it follows 
that / has a generalized convolution representation and (|9) holds. Conversely, assume 
that / has a generalized convolution representation and that (|9) holds. From (|9]l it 
follows that the Markov-parameters of yz and / coincide, i.e. S]^{sq) = SQ{sq) and 

S^j^ {qosq) = S'^j{qosq) for all s E Q* , q,qo E Q, j ^ I , . . . ,m. Since both yz and / admit 
a generalized convolution representation, by Remark|2]they are equal. The latter means 
that Lisa realization of /. □ 

Proof of Theorem^ It is enough to show that for any family of n x n matrices /vy, 
q E Q and any matrix G E M"^' for some / > the following holds. Define the matrix 



for e N. That is, is the span of the column vectors 

of FyG, V E 2^*^^^ Here we applied Notation |2] to Fg,q E Q to obtain the matrices 
Fv,v E Q*. Define the subspace as the space spanned by the column vectors of the 
matrices FyG, v E Q*. If we can show that lmMn-\ — ^ , then the statement of the 
theorem follows easily. 

Indeed, it is easy to see that the linear span of all reachable states of E equals 
J' , if we set Fq — Ag, q E Q and G = B. Moreover, in this case ^„_i — ^{T,). Hence, 
rank ^(E) = n is equivalent to ^ = M", which in turn is equivalent to span-reachability 
of E. Similarly, if we set F^ = and G = C^, then = ^„_i and is the 

orthogonal complement of flvee* kerCA,,. From |30| it follows that E is observable if 
and only if flreg* kerCA,. = {0}, which is equivalent to Im^„_i = = K". The latter 
is equivalent to rank i^(E) = n. 

We proceed to show — lm^„. The proof is the same as the one of an analogous 
statement for rational representations or state-affine systems 1, 18. .29,1 . We repeat it for 
the sake of completeness. It is easy to see that Im.^^. C for all ^ G N and ImMk C 
Im^<.+1 . By a dimensionality argument it follows that there exist < ^* < « — 1, such 
that Im^ii^ = lmMk,+i- From this, by noticing that Im^<.+i — ImG + ILqeQ ^^^q^k^ it 
follows that y — ImMk,. Since ImMk^ C Im^„_i, we then obtain that lmMn-\ — J" ■ 

□ 
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